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Abstract 

We consider a general formalism for treating a Hamiltonian (canonical) field theory with a 
spatial boundary. In this formalism essentially all functionals are differentiable from the very 
beginning and hence no improvement terms are needed. We introduce a new Poisson bracket 
which differs from the usual "bulk" Poisson bracket with a boundary term and show that the 
Jacobi identity is satisfied. The result is geometrized on an abstract world volume manifold. 
The method is suitable for studying systems with a spatial edge like the ones often considered 
in Chern-Simons theory and General Relativity. Finally, we discuss how the boundary terms 
may be related to the time ordering when quantizing. 
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1 Introduction 



Seen in the light of the renewed interest for theories where the edge plays a prominent role, cf. 
Maldacena's Conjecture Q, Carlip's and Strominger's different approaches for the microscopic 
counting of states on the (inner or outer) edge of the world |2|, 't Hooft's and Susskind's principle 
of holography |^], but also Chern-Simon theories and General Relativity ||5|, |^ in general, there 
is strikingly few papers devoted to develope the general formalism for Hamiltonian (canonical) 
field theory in the presence of a spatial boundary. Here we are thinking of the fact that the usual 
equal-time Poisson bracket 

fails to satisfy the Jacobi identity 

^ {F{t),{Git),Hm = 0, (1.2) 

cycl. F,G,H 

when space S has a boundary dT, ^ 0, at least if we apply the usual Euler-Lagrange formula for 



the functional derivatives in (1.1). (We shall show below how to ensure the differentiability of the 
functionals by using the notion of higher functional derivatives, so that the above violation is a 
fully legitimate problem to raise.) The failure of the Jacobi identity can be seen even in the most 
simple toy examples which have a non-trivial boundary. An equivalent manifestation of this fact is 
that functional derivatives cease to commute when a spatial boundary is present |0. 

The most common example of the above phenomenon is the usual d-dimensional flat space S = IR"^. 
Here the spatial inflnity \x\ = oo constitute a boundary for the space. This statement can be made 
precise by a so-called one-point compactification. 

An often used cure is to impose conditions on the dynamical fields (l)^{x,t) at the boundary which 
are consistent with the time-evolution. However, that approach may exclude interesting topological 
questions, such as solitonic field configurations. Our main goal in this paper is to see how far we 
can get without imposing boundary conditions. 

On the other hand, to calm the reader who perhaps finds these facts strange, let us mention that if 
"there is no boundary" (for instance, think of a torus, or equivalently periodic boundary conditions, 
or even vanishing boundary conditions), integrations by part inside the spatial integral does not 
produce boundary contributions, and the Jacobi identity for the above "bulk" Poisson bracket can 
be demonstrated after some straightforward manipulations. 



The paper is organized as follows: In next Subsection 1.2, we present a new Poisson bracket. 
Thereafter, we give some further introductional remarks about differentiability and improvement 
terms. In Section |2| we give a manifest formulation of the new Poisson bracket, discuss the higher 
Euler-Lagrange derivatives and develop a generator method (Section ^) which in particular is 
suitable of handling the arithmetic manipulations involved in the proof of the Jacobi identity. The 
proof itself is postponed to an Appendix ^. After that we turn to the questions that naturally 
arise with the existence of the new Poisson Bracket. Can it be given a geometrically covariant form 
(Section ^)? How does the boundary affect the Hamiltonian dynamics (Section ^)? To answer the 
last question, we have included a technical Section ^ to explain some supplementary formalism. 
Finally, we discuss the role of time order in connection with boundary terms (Section |^). 
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1.1 Notation 



The (f) {x,t), A = 1, . . . ,2N, denote the (bosonic) coordinate and momenta fields of the phase 
space. (Generahzations to fermionic variables are straightforward.) The non-degenerate symplectic 
structure is for simplicity taken to be ultra-local and constant. This vast simplification is the 
most interesting case for applications and already contain as we shall see an interesting structure 
worth analyzing by itself. Needless to say that a general field transformation <j) ^ cp' will violate 
such assumption. A manifestly covariant formalism under general field transformation is out of the 
scope of the present work. Also we assume, to avoid technicalities, that the d-dimensional space S 
can be covered by a single coordinate patch with a flat measure. (We shall relax these assumptions 
in Section]^.) Furthermore, in agreement with the spirit of the Hamiltonian (canonical) formalism, 
we shall assume that the functionals F{t) do not contain time derivatives {dt)^ (j)^{x^t) of the 
dynamical field variables cj)^{x,t). So if one for instance is interested in a (total) time derivative 
G{t) = d^F{t), where F{t) contains no time derivatives, one should study F{t) instead of G(t), etc. 
Finally, we assume that there is no temporal boundaries. This being said, time t plays no active 
role, and we can suppress the time variable t in what follows. 



1.2 New Poisson Bracket 



As mentioned in the introduction the bulk Poisson bracket (IT) does not satisfy the Jacobi identity. 
(A purist would perhaps then claim that ( |1.1[ ) does not qualify for being called a Poisson bracket 
at all! However, we shall continue to call it a Poisson bracket.) Knowing that the failure of the 
Poisson bracket ( |1 . 1| ) is at most a total derivative term, it is natural to speculate whether one can 
add a boundary contribution B{F, G) to this bulk Poisson bracket. 



= {F,G}(o)+i?(F,G) 



(1-3) 



so that the Jacobi identity is satisfied even in the presence of a boundary. In fact, this is so. We 
find that the following boundary term 



B{F,G) 



E 



(fx 



6F 



-UJ 



.AB 



6G 



(F^G) 



_6^M''){x) 6(p^{x)_ 
does the job. We have employed a multi-index notation: For instance the index 
k = iku...,kd)elN^\{{0,...,0)} , Wo = {0,1,2,...}^ 
runs over the d-dimensional non-negative integers (except the origo), and 



_d_ 

dx^ 



(L4) 



(L5) 



(1-6) 



(The main features of the construction are already present in the dimension d = 1 case. A first- 
time reader will not miss the essential points by treating the multi-index k as an integer, i.e. letting 
d = 1.) More importantly, in a perhaps conceptionally dangerous - but in practice convenient - 
notation, the 

denote the higher functional derivatives of F of order k. Here 



2 



is the usual functional derivative. In general, the higher functional derivatives are required to 
satisfy 



5F = (fx 



6F 



(1.9) 



for arbitrary infinitesimal variations (f)^{x) 
are not restricted at the boundary. 



(x) + 5(j)'^{x). In particular, the variations 5(j)^{x) 



A quick estimate shows that if the entries contain spatial derivatives di(j)'^{x) of the dynamical field 
variables cj)^{x) to order A^, then the full Poisson Bracket (|1.3| ) contains spatial derivatives up to 
order 3A^. Hence the algebra ^Ar=o of functionals with no spatial derivatives closes on the new 
Poisson bracket. However, physical interesting theories usually have functionals with up to = 1 
spatial derivative, i.e. they belong to the class An=i- This class An=i of functionals does not close 
on the new nor on the bulk Poisson bracket, and this is the main reason why we are force into 



summing over the index lattice (1.5) 



The idea of adding surface contribution is far from new. In a seminal work Regge and Teitelboim [|[ 
emphasized the importancy of having a boundary term in the action of canonical general relativity. 
However, they did strangely enough not add surface contributions to the Poisson bracket. Lewis, 
Marsden, Mongomery and Ratiu |^] considered a truncated version of ( |l.4D where only the terms 
with = 1 are present. Because they didn't add the |A;| > 1 terms, they needed to impose addi- 
tional boundary conditions to ensure the Jacobi identity. The first successful attempt to remedy 
this was conducted in 1993 by Soloviev |^. In our notation, his solution formula(3.4)] reads 



{F,G} 



E 

kl=0 



<fx d^+' 



6F 



-to 



AB 



6G 



6(l)^(''){x) 6(l)^i^){x) 



(1.10) 



It is easy see that his bracket is different from our solution (|1.4|). It would be interesting to know 



whether his bracket supports a manifest formulation (see Subsection 2.2 below), or more generally, 
if it is independent of the representative for the higher functional derivatives. After the first preprint 
of this paper appeared, Soloviev has made a comparison |1C] of the two solutions (|1.4|) and ( 1.10 ). 



We shall in this paper concentrate fully on the solution (1.4). 



1.3 Review of Differentiability and Improvement Terms 

The classical point of view H] on the problem with the Jacobi identity has been to view this as not 
so much a problem of the Poisson bracket itself, but rather that functional derivatives in general in 
the case of a non-trivial boundary dT, ^ are ill-defined when the functional, say F, depends on 
the spatial derivatives di(f>'^{x) of the dynamical field variables (l)'^{x). In this case there does not 
always exist functions /yi(x), such that the change in the functional F is fully described by 

= j d'^x fA{x) 6(l)^{x) , (1.11) 

for an arbitrary infinitesimal variation 5(j)^{x). (Of course in the affirmative case, we usually call 
fA{x) the functional derivatives of F.) 

Example. Consider an interval S = [a, b] and the toy functional 

F = tdxd^m = dm\i^a ■ (1-12) 
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The variation 5F can he identically rewritten as 



SF = j^{6j:{x,h)-5^{x,a))d54>{x) (1.13) 



To bring 6F of the form (1.11), one is tempted to do an integration by part. But this does not help 
us, partly because of the total derivative term does not vanish on the boundary. (Remember that 
we do not want to impose boundary conditions on the field (p, cf. the Introduction.) In fact, it is 
not apriori clear what should be meant (viewed separately) by any of the 2x2 terms arising from 
such an integration by part. So the traditional functional derivative is ill-defined. These difficulties 
should be compared with the ease that the same variation 6F is described by the higher functional 
derivatives ( [i..^ 

Si ■ (1.14) 



(5(/)W(a;) 



This is the general picture: A traditional functional derivative is often ill-defined or have a very sin- 
gular behavior at the boundary. (And at this point we haven't even touched the problems of building 



up the Poisson bracket (IT) out of two highly singular functional derivatives, i.e. multiplying two 
delta distributions together, both of which typically have support on the boundary. We should also 
mention that authors for such reasons often additionally require the functional derivatives f^ in 



Eq. ( 1.11| ) to be continuous functions.) In any case, this makes the traditional definition ( 1.11 ) not 



very suitable for systems with a boundary. 

Let us mention an important algebra of functionals, that are differentiable in this traditional 
sense (1.11), namely those functionals that do not depend on the spatial derivatives di4>^{x) of the 



dynamical field variables (j)'^{x). They form an algebra under the Poisson bracket. (In this algebra 
^0 the bulk Poisson bracket ( |1 . 1[ ) and the full Poisson bracket ( |1.3D coincide.) 



The traditional cure in case of an ill-defined derivative, is to improve the functional F -Pimpr 
with a boundary term, a so-called improvement term, so that the derivative becomes smooth and 
well-defined. The drawback is of course that we are then studying a different functional than we 
originally started out with! Typically, one meets the following preparation of an observable in the 
literature: A function f{x)=f{d^(f>{x),x) is smeared with a "test function" r]{x) into a functional 
of the type F[r/] = d'^x f{x) r]{x). It is then improved F[r]] ^ i^imprM ^ -^o so that it belongs 
to the above mentioned class Aq by recasting all the spatial derivatives to hit the test function. 

We will bypass all this, i.e. the bottle-neck of ( |1.11| ), by using more functions (the higher functional 
derivatives) in the description of an arbitrary variation SF. The format of (^]^) is so broad, that it 
essentially covers all interesting functionals, which do not contain time-derivatives, cf. the discussion 
in Subsection |l.l[ (One could of course give meaning to differentiation of a functional with temporal 
derivatives simply by brute force extending the multi-index k in (|1.9| ) from d dimensions to d -|- 1 
dimensions. Although relevant for so-called covariant formulations (covariant in the sense that time 
and space are treated on equal footing), this is not in the line of the Hamiltonian theories, and 
hence not something we will pursue in this paper.) 



2 General Formalism 
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2.1 Partial Derivatives of a Functional 



Let us describe the higher partial derivatives of a functional F. (They are not to be confused with 
the usual higher partial derivatives of a function, although they are related.) In fact, they are 
objects, given the suggestive notation 

dF 

d^Mk){x) ' ^^"^^ 

that satisfies 

dF 



6F 



r °° BF 

/ d'^^ y ^ , d^5<i)^{x) , (2.2) 



for arbitrary variations Scf) {x). If this notation (2.1) in the future leads to ambiguities, we will 



specify whether we mean partial differentiation wrt. a function or a functional. Usually the context 
will exclude one of the possibilities. In fact, in this article we will often use the notation 

^..^)/(-) = ^0^) ■ (2.3) 
for the usual (higher) partial derivative for a function f{x) = f {d^ (/}{x) , x) . 
2.2 Manifest Formulation 

The set of the higher functional (and the partial) derivatives may not be uniquely defined, so one 



may worry that the full Poisson bracket given by (O) and (L4) depends on the choice of the 
representatives for the higher functional derivatives. The answer is that it is independent. This 
follows from the manifest formulation given below. 

We begin by giving a more useful definition of the (usual) functional derivatives (^]^) than the 
traditional definition, cf. Eq. ( |1 . 11 ) . The differential 5F = 5F[(f), 6(j)] of F is assumed to split 



5F = dF + dF (2.4) 

into a bulk integral 

r r)F 

dF = dF[4>,6<f>] = / d''x—^6^^{x) (2.5) 

and a boundary integral 

d „ 

OF = dF[(l),6^] = y d'^x di F\5(t)^{x) , (2.6) 

where F\ in general can be differential operators acting on 5(j)^. That the integral (^]^) is a 
boundary integral is justified by the divergence theorem. If we furthermore require the (usual) 
functional derivative to be a contineous function, this function together with the above spilt (^^) are 
uniquely defined. (We stress that the uniqueness of the (usual) functional derivative is jeopardized 
if for instance it contains a delta distribution with support on the boundary. This is ruled out by 
requiring continuity. See also the uniqueness discussion in the next Section.) 

The bulk Poisson bracket {•, •}(o) is now well-defined by the Eq. ( |1.1D , and the full Poisson bracket 
Eq. (|l.3| ) differs from this by a boundary term 

B{F,G) = 5fU,<50={0,G}(o)1 - {F^G). (2.7) 
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At this point the reader can merely take 5(/) = {0, G}(o) as being a convenient short hand notation 
for 



There is two obvious ansatzes for the differential operator F^. Either with the spatial derivatives 
ordered to the right of some coefficient functions (also called the normal order of x and d). Or 
vice- versa. In the former case 

we call the coefficient functions for the higher partial vector derivatives. The name "vector" refers 
to the i-index. In the latter case 

n - Zfj^ . (-0) 

the coefficient functions are called the higher functional vector derivatives. 



2.3 Uniqueness of the Higher Derivatives 

Up to now, we have only characterized the higher functional (partial, vector) derivatives of a 
functional F in a descriptive manner. The question of existence yields rather mild conditions, that 
we shall not be concerned about. The issue of uniqueness is a much more interesting question. The 
ambiguity in the choice is most clearly displayed via the higher partial vector derivatives: 

Uniqueness Theorem. Assume that there is given a sequence of contineous functions f^ and 
f{k)A ^^'^^ vanish identically with the exception of a finite number and such that for an arbitrary 
variation 6(f) 

„ d „ oo 

= d-'x fA{x) 5cp^{x) + d^xd^Y. /^(fc)(x) d''5cl)^{x) . (2.11) 

Then fA = in the entire space S and the (higher) partial vector derivatives f\(^f^-^ are tangential 
to the boundary dT,. In detail, 

VxeaS : /A(fc)(x)^ (/A(fe)(x),...,/i(fc)(x)) eT,(aS) . (2.12) 

That the first function /a = vanishes is just a restatement of the uniqueness of the (usual) 
functional derivative, while the boundary condition refiects that the higher partial vector derivatives 
may be modified with a vector field that doesn't locally carry a boundary fiux. 



2.4 Relations among the Different Kinds of Higher Derivatives 

In order to get an idea of how ambiguous the other variational discriptions are, let us give some 
maps between the mentioned choices of the higher functional (partial, vector) derivatives. We start 
with a bijective correspondance between the two scalar definitions ( |1.9D and (|2.2| ) of the higher 
derivatives. If there is given a sequence of the higher partial derivatives satisfying the definition 
(O), then 






(2.13) 
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satisfies the definition (p. .91). On the other hand, starting from the point of view of the higher 
functional derivatives, we get a solution to the higher partial derivatives by 



dF 



E 

m>k 



The variational descriptions ( |1.9| ) and 



m 
k 



d 



,m—k 



6F 



dF 



coincides because of the x-pointwise identity 
6F 



(2.14) 



A:=0 



-50^ (x) 



(2.15) 



For a proof, see the equation ( [A.2| ) in the Appendix. Similarly, one may transform bijectively 
between the two vectorial definitions by use of the corresponding relations 



6'F 



(5(/.^W(x) 



E 



m>k 



m 

k 



i-d) 



m—k 



d'F 



E 



m>k 



m 
k 



m~k 



6'F 



(5(?i^M(x) ■ 
(2.16) 

Not surprisingly, as the vectorial definitions (|2.g| ) and ( p. 10 ) carry the most indices, they have the 
greatest flexibility in representing a solutions. We may convert from the higher vector to the higher 
scalar derivatives via the formulas (for k ^ 0): 



dF 



a(/)^(^)(3 



E 

ielik) 



:+E3' 



5F 



5^F 



d(l)Mk-e,)(x) fr{ 9(/>^W(x) ' 6cl)Mk){x) 



Sk) s^^^'-^^Kx) 



(2.17) 



Here = (0, . . . , 0, 1, 0, . . . , 0) is the z'th unit vector in the index lattice and I{k) = {i = 1, . . . ,d \ ki ^ 0}. 
Going from the higher scalar to the higher vector derivatives is not unique. A natural choice is for 
the functional derivatives are 



5^F 



1 



5F 



50^W(x) \I{k + ei)\ <5(/>^('=+eO(x) 



(2.18) 



We leave out the corresponding relation between the partial derivatives to carry on with our main 
application: Local field theories. 



2.5 Local Field Theory 

Let us restrict ourselves to local field theories, i.e. all functionals can be expressed as an integral 
F = l^d'^xfix), fix) = f[d'cP{x),x) , (2.19) 

where k runs over a finite subset of INq. Note that we have allowed for explicit x-dependence in 
/. It essentially costs no extra work, and it becomes important later on. We shall postpone the 
analysis of functionals that depends on external space-points to a later section (Section ^), partly 
because it would be notational inconvenient to address those now. 



As mentioned before the (higher) partial derivatives of F need not be unique. Our strategy will 
be to "lower" the definition from the level of integrals to the level of integral kernels. In other 
words, this means that if a functional has different integral kernel representations, this may lead 
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to different definitions of the higher derivatives. For instance, at the case at hand, i.e. of a local 
functional ( p.l9 ) with a given integral kernel /, there is one natural candidate 



dF 



PA{k)f{x) , 



(2.20) 



that fits the relation (2^). So strictly speaking, the (higher) partial derivatives are really (higher) 
partial derivatives of the kernels, although we will not indicate this explicitly in the notation. In 
the same fashion the distinguished candidate for the functional derivatives becomes the higher 
Euler-Lagrange derivatives: 



5F 



EA(k)f{x) = E ( 

m>k \ 



m 

k 



i-d) 



m—k ] 



A{m) 



fix) . 



(2.21) 



For a mathematical textbook on higher Euler-Lagrange derivative, see for instance Olver |11, p.365- 
367]. Note that the m-summation in ( p. 21 ) terminates after finite many terms in case of a local 



functional, so that we do not have to worry about convergence of the sum. Let us simply use (2.20) 
(and ( |2.2lD ) as the working definitions for the local functionals. After all, our primary goal is to 
prove the Jacobi identity for the full Poisson Bracket, and this does not depend on the choice of 
the representatives for the higher derivatives. It is evident from the natural solution ( 2.20]) (and 
( p. 21 )) of the (higher) derivatives, that in practice all the local functionals that one encounters in 
physics are differentiable. 

2.6 Ultra-local Poisson Bracket 

Having restricted ourselves to the ultra-local case, let us for each pair of local functional 



F = I (fx fix) and G = f (fx g{a 
define a x-pointwise version 

{f.9}{x) = {/,<?}(o)(x) + i?(/,5)(x) 



of the Poisson bracket 



{F,G} = J^d''x{f,g}ix) = {F,G}(o)+i3(F,G) . 



(2.22) 
(2.23) 
(2.24) 



Namely define 



{f,9}{o){x) = EA(^Q^f{x) uj"^^ EB{o)g{x) 
B{f,g){x) ^ 5] [i?^(fe)/(x) cu^^ EB(0)5(a 

This means that the full x-pointwise Poisson bracket reads 



if^g) 



(2.25) 



{f,g}{x) = '9'[^A(fc)/(x)^^^i^B(o)5(x)] --{/,5}(0)(a;) - if o) 

k=0 

= E PAik)fix)uj^'' d''Es(o)9ix) -::{!, g}(o)ix) - (/ - 5) • 



(2.26) 



k=0 



The last equality in ( 2.26 ) follows from equation ( |2.15| ). We can now conduct our analysis x- 
pointwisely. In the next Section H we shall suppress the space point x € S. 
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3 Generator Methods 



3.1 Heisenberg Algebra 

Due to the quite heavy combinatorics involved in the proof of the Jacobi identity, it is useful to map 
the above problem into a simpler and in fact well-studied object, namely the Heisenberg algebra. 
Although the actual proof is presented in the Appendix, we find the central idea, while perhaps 
not entirely original, is quite important, so we will present it here. For a recent exposition of Fock 
space methods for variational systems, see also |T^ ]. 

Let us study the interplay between the partial and the spatial derivatives. The higher partial 
derivatives PA{k) commute among each other, but they do not commute with the (total) spatial 
derivatives 

= ^^(^•+^')p^(fc)+5rp""* , (3.1) 

where = (0, . . . , 0, 1, 0, . . . , 0) is the i'th unit vector in the index lattice. More precisely, we have 

min(fc,n) / \ 

PAik) = E (Z) ^""""PAik-m) . (3.2) 

m=0 \ I 

The main idea is to simulate this complicated disentanglement formula with the help of a Heisenberg 
algebra. Let us introduce abstract (bosonic) algebra elements Y\ that obey the following Heisenberg 
algebra commutator relations 

[yi,a,] = 5;. , [yi,y^] = o, [5„9,] = o. (3.3) 



The third equation is not a definition, but is a well-known consequence of (3.1). It is a remarkable 
fact that we can mimick the non-commuting behavior of formula (|3.2D by formally writing the 
higher partial derivatives product 

PAik) - p7-^_ (3.4) 

of the Y\ algebra elements and what is basically reduced to be a passive spectator in what follows, 
namely 

We take Pa to commute with everything: 

\PA.yh\ = 0, [PA,^^] = 0, [Pa,Pb] = 0. (3.6) 
Above we have adapted the following multi-index conventions: 

Ya = (Yl)"' ■ ■ ■ {Yi'f\ k\ = h\---kd\. (3.7) 



That the Heisenberg algebra ( |3.3| ) with the formal assignment ( p.4| ) really reproduces (3.2) is proven 



in the Appendix see equation ( |A.1| ) . The proof becomes very simple once we adapt the generator 
techniques of the next section. 
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3.2 Generator Methods 



As a second computational improvement, it is useful to hide the integer indices inside generating 
functions which depend on continuous parameters qt, i.e. we shall sum up in generalized Fourier 
series. We implement this program for the (higher) partial and the (higher) Euler-Lagrange deriva- 
tives, respectively, as follows 



PAiq) ^ Paw*?' = Pa e"^-^ 



EA{q) = EA(k)Q = ■.exp[{q-d)YA]PA- = exp 



-d— 

dq 



PAiq) ■ (3.8) 



Prom now on we will implicitly imply summation over repeated multi-indices k € INq. In fact, we 
may view the multi-indices sums as running over the entire d-dimensional integer lattice k € 
by simply declaring that objects like 



Yi, , PA(k) , EA,k) , (3.9) 



are zero if k is outside the original non-negative d-dimensional quadrant INq. The next-to-last 
equality in ( |3.8| ) follows from the mere definition of the (higher) Euler-Lagrange derivatives (2.21), 
once we have declared the following normal ordering prescription: 

■.YXdr- = -.d.YX: = d,YX. (3.10) 

3.3 Fourier Transform 

As a third computational improvement, let us Fourier transform the variables qi to variables y*. 
PA{y) ^ I d^q e-iyPAiq) = Pa 5(1a - v) , 

EAiy) ^ J d\ e-'^ySAiq) = e-^^P^(y) = Pa 5(5^^) e"^^ . (3.11) 

With the theory developed so far, we have reached the second objective of this paper (the first 
objective being to give the form of the full Poisson bracket ( |l.3D ), namely achieved a formalism that 
are capable of giving a short proof of the Jacobi identity. For the proof itself, see the Appendix |B|. 
Note that because the Poisson bracket is independent of the actual choice of representatives for 
the higher functional (partial) derivatives, it implies no limitation that we use the natural choice 
( p. 21 ) and ( 2.20| ). We will now turn to the question of geometrizing the Poisson bracket. 



4 Abstract Manifolds 

In this section we formulate the results obtained so far in a geometrically covariant manner in- 
dependent of the choice of coordinates. More precisely, the construction is generalized from a 
d-dimensional subset S C M"^ (where the space and the chart are identified) to an abstract d- 
dimensional manifold S with spatial covariant derivatives Di and a d-dimensional volume density 
p. In other words, the spatial derivatives di are replaced by covariant derivatives (let us indicate this 
in an oversimplified way as Di = di + Pj), and the trivial measure d'^x is replaced by pd'^x. We will 
assume that Di = Di(x) and p = p{x) do not depend on the dynamical fields 0'^(x, t) nor on time t. 
We do not assume that the volume density is covariantly preserved, i.e. that Dip = {di — ^^k)p = 0. 
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In passing from derivatives di to covariant derivatives Di , we face the main complication compared 
to the flat case. In general, the spatial covariant derivatives do not commute, when the curvature 
is non-vanishing. We have 

[Di,Dj] = [£)^^xplicit^ ^explicit J ^ ^4 -^^ 

where the (total) covariant derivative Di is given as 



D, 



'^A{K) + Di 



(4.2) 



The index structure of the first term in (|4.2| ) will be explained below. In the general case of non- 
vanishing curvature, one can proceed by declaring that the functionals of the theory depend on 
ordered tuples of covariant derivatives (jy'^^^^x) = D^(j)'^{x) of the dynamical fields (j)"^{x), rather 
than only unordered sets of derivatives (/i^^'^) (x) . An ordered tuple K is of the form 



K 



{ki,...,k\K\) G {!,..., d} X ... X {!,..., d} , 



(4.3) 



where d is the space dimension, i.e. the dimension of S. We have given a resume of the calculus of 
ordered tuples in the Appendix ^ All formulas carries over to the curved case in essentially the 
same format. However, there are some notable differences that we now stress. The description of 
the higher derivatives (|2.2D, (1.9), ( 2.13| ) and ( |2.14| ) are replaced with 



6F 



OF 



SF 



dF 

a0^(^)(x) 



p{x)d'^X ^ r.,AlK)( 

f p{x)d^x £ Z)^ ^ 



S(j)A{K)(^X 

dF 



D^6(t>^{x) 
6^^{x) 



M>K 



dct)^(^''){x) ' 

5F 



(4.4) 



For details concerning the notation in ( [1.4[) , see definitions in Appendix p.l| . The (higher) functional 
and partial derivatives inherits tensor properties, if F is covariant. So the formulas are covariant. 
The bulk and the boundary term of the Poisson bracket reads 



{i^,G}(o) 



p{x)d'^x 



6F 



-UJ 



AB 



6G 



6^^{x) 5(1)^ {x) 



B{F,G) = ^ ^ />(x)/x L*^^ 



6F 



.AB 



6G 



(F^G) 



(4.5) 



Note the apparent asymmetry between the two last formulas in ( [4.4j ) with a transposition t of the 
order of covariant derivatives in the third equation. As a rule of thumb one may say that the spatial 
derivatives {—d)^ in the flat metric formulation becomes {—D)^ in the covariant formulation. This 
generic feature carries over to the generator formalism (p.S|): 



PA{q) = Pa{K)(1 
EA{q) = EA{K)q 



K 



{-DY 



exp 



exp 

d_ 
dq 



D 



d_ 
dq 

PA{q) . 



EAiq) , 



(4.6) 



Let us note that the g^'s (besides commuting with everything else) does not commute among 
themselves. More precisely, they are freely generated. This is necessary in order not to loose 
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information about the operator ordering when passing to the generating functions. Also the Fourier 
transform ^ yi, cf. Section |3.3|, can be given sense in the non-commutative case. 



The replacement of the disentanglement formula (3.2) becomes 



M<K,M^N 

Pmk) D"" = J2 D^^^'PAi-M+K) . (4.7) 

We are able to define contravariant elements ^^(A')! such that 

Pmk) ^ Pa^- (4.8) 
The Heisenberg algebra ( |3.3| ) is traded for 

[y4(x),A] = \K\ YA^-i+K) ■ (4.9) 

Remarkably, even in this non-commutative case, the exponentiated version can be recasted into the 
following simple form 

We shall have more to say about this construction at the end of Appendix The main point is 
that the proof of the covariant Jacobi identity can be demonstrated in almost exactly the same 
way as in the flat metric case, cf. Appendix 

Note that for non-zero curvature [Di, Dj](f)'^{x,t) ^ 0, the actual field value 4>^{x,t) is apparently 
not well-defined, i.e. if one tries to sum up the change in (l)^{x,t) along a closed loop, one obtain 
a non-zero result. This is worse than a global obstruction. Perhaps it should be called a local 
obstruction. (A similar situation occures, say, in bosonic string theory with the Polyakov action 
when the worldsheet metric has a non-zero curvature. This also leads to problems in locally 
assigning values to the target space fields.) The problem seems less formidable in the context of 
the Feynman path integral, where we only assign field values along one path at the time. But it 
is a genuire challenge for the operator formalism. One way of making sense out of this would be 
to declare the decendent fields (p^^^\x,t) to be independent fields living in a non-commutative 
jet-bundle. In any case, we feel that it would be too hasty a priori to draw conclusions in general, 
and we leave it to the future to appropriately implement non-zero curvature in specific physical 
theories. 



5 Supplementary Formalism 

Until now, we have only discussed functionals F with no external space dependence, i.e. all space- 
variables are integrated out. However for physical applications, we would like to conduct manipu- 
lations directly on the integral kernels rather than the integrals. For instance, to give sense to the 
fundamental equal-time relations 

{<P^{x,t),cl>^{y,t)} = a;^^(x,y) = u""^ 5j,{x,y) . (5.1) 

The plan for the rest of this article are 

• to treat the Dirac delta distributions (and the derivatives thereof) in the presence of a bound- 
ary. Distributions is a vast subject in their own right, and we will here only give a heuristic 
treatment. 
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• to extend the definition of the higher functional derivatives to more general types of local 
functional. 

• to analyse the implications for the Hamiltonian dynamics. 

While parts of this section is standard material, it is reviewed for continuity and to fix notation. 

5.1 Embedded Approach 

Having a geometrically covariant formulation at our disposal enables us to reduce the discussion 
to a single chart. We can slice up space in smaller regions; thereby producing unphysical double- 
sided boundaries (unphysical domain- walls) , and we can hence consider space within such a smaller 
region S covered by a single chart. The local geometric data about the physical space S is stored 
in the volume density p and the covariant derivatives -Dj. 

Furthermore, we will assume that S takes place inside a bounded region of the chart IR!^, i.e. that 
it can be placed inside a large ball in IR'^. Note that we are not placing any restriction on the 
distances in the physical space S; only on the distances in the chart. Or perhaps we should say: 
in the choice of the chart. For instance, if space E = IR^ is the ordenary flat space, one should 
map flat space into a bounded region S of the chart IRf" using a non-trivial p and Di. In this case 
the spatial infinity is truly the boundary of the region S. The perspective will be that of a typical 
Penrose diagram: "There is always room for something beyond spatial infinity." The motivation 
for the above assumption is deeply founded in the theory of distributions, cf. below. 

This being said, we will adapt the usual practice of identifying the space S C IR!^ with a region of 
the chart. 

5.2 Reguleirized Chciracteristic Functions 

Let us consider the characteristic function 

^ f 1 if iceS 
1 otherwise , 

for the space S C ]R'^ as a limit of a smooth function Xe(x), where < Xf^[x)^lT,{x) for e — > 0+. We 
regard Xg = X^{x) as independent of the dynamical fields 4>^{x) and as a scalar under coordinate 
changes x x' . The actual implementation of the regularization X^[x) should not matter, so one 
might as well choose a convenient form. One could for instance do as follows. Let dY.{x) denote the 
signed distance from x to the boundary C JR*^ as measured in the chart IR'^. The signed distance 
di:{x) > is positive if a; G S° belongs to the interior and it is negative dY,{x) < if x G (iR'^\S)°) 
belongs to the exterior. Then we could implement Xe(x) as 

(5.3) 

(This looks horrible in other coordinates, so a geometrically minded reader might prefer to sub- 
stitute the chart IR'^ with an abstract unphysical embedding manifold. We shall not explore this 
point of view further in this paper.) 



X,{x) 



1 -I- exp 



dT.{x) 
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Next, we extend p and Di smoothly (and arbitrarily) to the unphysical sector lR'^\Ti. It may happen 
that p or Di themselves are singular at the boundary 9S. In that case one should consider smooth 
regularized functions p^ or DI, that in the limit e — > reproduces p and D^. Then an integral over 
S should be though of limit 



p{x) (fx fix) = lim / X,{x)p,{x) (fx f,{x) . (5.4) 

where the integrand f{x) also should be smoothly (and arbitrarily) extended to the unphysical 
sector 1R'^\T,. (In case of more than one region, X^{x) should be a differentiable partition of the 
unity.) From now on we will not write nor question the limit e — > 0, but merely take for granted 
that this e-prescription makes sense. 

5.3 Dirac Delta Distributions 

Throughout the paper, the Dirac delta distribution 6{x — z) refers to the full iR'^-chart, while (5s (x, z) 
refers to the physical space S. The physical Dirac delta distribution (5s(x,z) is characterized by 
the property: 

Vr/ : / p{x) (fx p{z) (fz rj{x,z) 5j]{x,z) = I p{x) (fx rj{x,x) . (5.5) 

We can realize the physical Dirac delta distribution 5s{x, z) in terms of the unphysical Dirac delta 
distribution 6{x — z) as 

The main idea behind demanding that S should occupy a bounded region of the chart M'^, is that we 
can then perform formal integrations by part on the unphysical Dirac delta distributions 6{x — z). 
This is so because we can consider all test functions as having a bounded support in the chart M'^. 
"Test functions" should here be read in the broad sense of the word that in particular includes 
functions f {(p^^^x) , x) of the dynamical fields (j)^^^\x). 

On the other hand, integration by part of the physical Dirac delta distribution 6j:,{x,z) will in 
general lead to boundary contributions at the physical boundary dT,. The detailed form can be 
inferred from the above relation ( |5.6| ). The benefit of this procedure, is that we do not have to 
postulate perculiar rules for the physical delta distribution. They may simply be derived from 

(El)- 

The above is the key observation in our analysis of distributions. Mathematicians have always 
(and presumably for good reasons) considered test functions in IR'^ as having compact support. We 
observe that if space S, which itself could be unbounded, fills a bounded region of the chart M"^, 
we can without touching the above principle, still probe boundary issues at the physical boundary 

as. 



From (5^) and (|5.q ) it also becomes clear that the study of a non-trivial boundary dH and the study 
of a non-trivial volume density p are intimately related. With the spatial integration interpreted 
as ( |5.4D , we may define the adjoint Dj of Di by formal integration by part in the chart M'^. It 
becomes 

Dj = - ^AXep(-) . (5.7) 

A^p 

(The arrow over Di indicates that the derivative Di acts all the way to the right.) 
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5.4 Space of Functionals 

Consider now a function depending on variables ^(i) , . . . , Z(^r) ; 

fiz^,y...,Z^r)) = , (5.8) 

where , . . . , K^^^ are multi-indices. For convenience we shall often use the compact notation 
z = (2;(i) , . . . , Z(^r) ) if there is several space points . We shall restrict ourselves to the space A of 
functionals F that can be expressed as a s-fold multiple integral over E 

Fiz) = I pd^z^^y.-pd^z^,) fiz) , (5.9) 
JSx...xS 

for some s G W. (It is implicitly understood that the 2;(j)'s which are integrated out on the right 
hand side, do not enter the argument list on the left hand side.) Furthermore, this space ^ is a 
(T-vector space. It is stabile under multiplication whenever defined. (Recall that the product of 
distributions need not be well-defined.) It is closed under integrating out external variables, or 
identifying external variables, say = Z(^jy We shall see below that it is also closed under the full 
Poisson bracket. 

5.5 Suitable Form of Functional 

Consider a local functional F{z) = J (i'^X(i) . . . d'^X(^r) fi^, z) with external dependence z. The typ- 
ical integrand /(x, z) consists of 

1. The Dirac delta distributions 5(x — y) and 5Y,ix,y)\ The regularized characteristic functions 

(x) in some power p E M. 

2. Smooth test functions with bounded support. This in particular includes smooth functions 
gi(f)^^\x), x) of the dynamical fields (p'^^^^x) and the smooth volume density pix). 

3. Derivatives Di acting on the various factors of the integrand mentioned under point 1 — 2. 
The above listed objects appears in two versions: 

A. an external type, if it depends on non-integrated external variables. 

B. an internal type, if it (at least partially) depends on integrated internal variables. 

An integral (kernel) is declared to be on suitable form if all internal derivatives (-B3) act on type 
B2 objects. In other words that the more singular type (Bl) objects are not hit by the the internal 
derivatives (S3). 

A functional (kernel) /(x, z) is not well-defined if one cannot obtain a suitable form by pure algebraic 
manipulations. In practice, this means 

• after formal internal integration by part, 

• after use of the Leibnitz rule and linarity. 

• after use of the identity igix) — giy)) Six—y) = , 

• and after use of the identity (-D^j.) + -D(j^) J 6ix—y) = 
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for the (partially) internal variables. 



Just from the freedom to perform a formal integration by part on the internal differentiation {B3) 
(or choosing not to do so, respectively), there is 2" ways of writing down an functional, where n 
is the number of internal differentiations (-B3). In practice, in all interesting functionals, every 
internal variable Xj appear at least once in the argument list of an internal singular object {i.e. of 
type Bl). As a consequence, in these cases, at most one of the above mentioned 2" choices leads 
to a suitable form. 



Needless to say that if one also integrate out the external z-variables F(z) without yielding enough 
room for the smearing type 2 objects, the result may not be well-defined. 



5.6 Higher Partial Derivatives 

Consider a functional F{z) € ^ in the algebra A. Assume that the functional (kernel) is of suitable 
form. Let us now define the higher partial derivatives as 



I pd^z^^) . . . pd'^z^,) ^^(^^ ho) PAiK)fi^) ■ (5-10) 

•^^^ i=l 

F{z) = f[d''cP{z),z) , (5.11) 

= E^E(x,Z(,))pi(«V(.). (5.12) 

We can formally extend the application range of the above equation ( |5.12| ) to include functionals 
F{z) as well, by implicitly assuming that the internal delta distributions automatically are placed 
inside the integration symbol. Then ( |5.12| ) becomes a convenient shorthand notation for ( ^.lOj ). 



a0^W(x) 

In case of a function 
this reduces to 



5.7 Higher Functional Derivatives 

The general definition ( [4. 4]) for a functional of suitable form yields 

i=l 

E E {-D^.)r'^'''Sj,ix,z^,^) P'll^.F{z) . (5.13) 



Ea(k) E [^^ Hi) ) ^ (^(1) ' • • • ' ' ^' ' 



A{My 

i=l M^K 



The above derivatives of a delta distribution may be resolved in two different ways: 

• By inner evaluation: The derivatives leave the delta distribution 5Y.{x,Z{i^) via the 

leg. If there is enough internal integrations inside the functional F, one may resolve the 
delta distributions by integration, thereby prolonging the derivative to an object inside the 
functional F. If all terms are to be resolved this way, this means that all the ^-variables have 
to be internal. 

• By outer evaluation: The derivatives leave the delta distribution 5y,{x, zi^i-^) via the x-leg. 
We await an external integration over the x-variable to evaluate the derivative of the delta 
distribution by formal integration by part. Let us stress the fact, that if one rely on the latter 
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method for the Sy:{x, term in ( 5.13| ) with the being an internal variable for the F(z) 
functional, then the x-integrated version is not on a suitable form as it stands. 

This distinction is important if one is to conduct further partial differentiations wrt. to the dynam- 
ical fields (j)^{x) on the functional (kernel). However, if no further differentiations are performed, 
the two methods yields the same result. 

As the most important example, we mention 

Note that the above definitions ( |5.12| ) and ( 5.13| ) guarantee the linearity and the Leibnitz' rule of 
the (higher) partial and functional derivatives. We also find that two (higher) partial derivatives 
commute. One may show in the case of a vanishing boundary dT, = 0, that the usual functional 
derivatives commute. In the case of a non-trivial boundary dT, ^ 0, the higher functional derivatives 
(as well as the usual functional derivatives) do not commute in general. 

5.8 Induced Functional Derivative 

Finally, one may define a induced functional derivative from the perspective of the embedding 
manifold, i.e. the chart M'^: 

for an arbitrary variation 6(f)^{x). This is so, because S C ]R'^ is bounded inside the chart M'^, so 
integration by part yields no boundary contributions at |a;| = oo. The induced functional derivative 
makes sense, eventhough there appears coinciding space points, because X^{x)p{x) does not depend 
on the dynamical fields (j)^{x). We can write it constructively as 

6[X,{x)p{x)F{z)] ' 
It is related to the higher functional derivatives via 

This induced functional derivative has the remarkable property of commuting with the spatial 
derivatives 

Ah,,) 5[X,{x)p{x)F{z)\ ^ 5[X,{x)p{x)Df'''h{z)] 

5(l>A{x) 6cP^{x) ■ ^ ^ 

This should be compared with the corresponding behaviour of the usual functional derivatives: 

' 6(f>^{x) 64>\x) ■ ^ ' 

The induced functional derivative satisfies the Leibnitz rule and it commutes with integrations: 



X,(x)p(x)/2pd'^^(i)F(2) 



, 5[X,{x)p{x)F{z)] L - - -J .,^n^ 
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5.9 Annihilation Principle 



As mentioned before, integration without smearing can produce ill-defined terms. However it can 
be very cumbersome to a priori discard all the bad terms of an expression. We shall therefore 
formally allow ill-defined terms to appear by giving a prescription, that consistently identify them 
and put them to zero. This is done by defining a little more restrictive version of the above so-called 
suitable form. Notation: For simplicity, we will assume from now on that the covariant derivatives 
commute, i.e. that the curvature vanishes. Assume further from now on that the volume density p 
is covariantly preserved. Dip = 0. A typical functional (kernel) consists of 

1. Dirac delta distributions 5{x — y). 

2. Regularized characteristic functions X^(x). (Integral powers X"(x), n > 2 should be considered 
clS 3j Tl- fold product of elementary X^(x).) 

3. Negative powers X~p{x), p > 0, of the regularized characteristic function. 

4. Smooth test functions with bounded support. This includes smooth functions g{(f>^^\x),x) 
of the dynamical fields (f)^^^\x) and the smooth volume density p{x). 

5. Derivatives Di acting on various factors of the integrand mentioned under point 1 — 4. 

We may assume by use of the Leibnitz rule and breaking the integral F{z) into several terms if 
necessary, that all derivatives {A5 + B5) only act on one elementary object under point 1, 2 and 4. 



(Here the letters A and B refers to the notation introduced in Subsection 5.5, while the numbers 
1 — 5 are those defined above in this subsection.) 

A functional F{z) of the above atomic type is declared to be identical zero if one cannot by algebraic 
means, cf. above, obtain a form where all internal derivatives (-B5) acts on type Bi objects. (Or 
in other words, the more singular type of objects Bl — B3 are not hit by the derivatives.) 



6 Hamiltonian Edge Dynamics 

In this section we discuss implications of the full Poisson bracket for the Hamiltonian dynamics. 
We first have to extend the definition ( |4.5| ) of the full Poisson bracket to more general functionals 
with external dependence. As a first principle for writing down the more general Poisson bracket, 
we shall demand that integrations f^.P^'^^ii) commute with the Poisson bracket {•,•}, that is 

{F{z),G{w)] = {J^pd^z(^,^ F{z),Giw)} . (6.1) 

This principle leads naturally (modulo the action of the annihilation principle) to what we shall 
call the solid Poisson bracket. We shall later see that it can be recasted into a so-called floating 
Poisson bracket, that (at a superficial level) takes different shape on different types of functionals. 
However, one may show by applying the annihilation principle that no actual differences take place. 

6.1 Solid Poisson Bracket 

Using the extrapolation of the formulas in the previous sections, the full Poisson bracket becomes 

{F{z),Giw)} = {Fiz),Giw)}^o)+B{F{z),G{w)) , (6.2) 
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where 



{F{z),G{w)\o) 



B{F{z),G{w)) 



p{x)(fx 



6F{z) ab5G{w) 
5(t)^{x) 



ExS 



p{x)d'^x p{y)d'^y 

AB ^G{w) 



- (F{z) ^ Giw)) . 

Alternatively, we may write the full Poisson bracket as 

dF{z) 



6(p^(x) 
5F{z) 



-LJ 



{F{z) ^ G{w)) 



dF{z) 



AB 



SGjw) 
6(1)^ (x) 



{F{z),G{w)} 



p{x)d'^x p{y)d'^y . ^ 



■UJ 



A{M)B{N) 



{x,y) 



dGjw) 
d(j)BW{y) 



(6.3) 



(6.4) 



where the symplectic kernel oj^^^^'^^^^^x ,y) reads 



UJ 



A{M)B{N) 



x,y) 



UJ 



.AB 



-D 



t ^^n 
(x)) y 



-D 



t \N 



Ssix,y) 



UJ 



AB 



+ ^--^(x)) ^(y)+ ^{x)y-^{y)) 



Sj:ix,y) . 



3.5) 



In the case where (at least) one of the M and are 0, the symplectic kernel can neatly be written 
as 



M 



V N 



UJ 



A(M)B{N) 



Note on the other hand, that the case M = N 



ix,y) = u;^^Df^)Z)(;)fe(x,y) . (6.6) 
is also the maximal case to make fully sense 



out of the expression D^^s^D^^6^{x,y) without employing the annihilation principle. Beyond that 
case, i.e. if Af 7^ A A 7^ 0, there is no escape ways left open for the X~ -^-function (via formal 
integrations by part of the derivatives). It is sandwiched inside the delta distribution between 
derivatives. It should be merged with a X^-function outside, whose mere existence on the other 
hand prohibes that a suitable (and hence a well-defined) form can be reached by means of integration 
by part. 



6.2 Hamilton Equations of Motion 

Consider the Hamilton equations of motion 

d 



dt 



-F z) 



{H,F{z)} + ^F{z) 



(6.7) 



Let the Poisson bracket be the solid Poisson bracket ( |6.3| ). And H = Jj, p{x)d'^x Ti.{x) be a local 
Hamiltonian. Then the Hamilton equations of motion for the fundamental fields read 



d_ 

di^ 



UJ 



AB 
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We also get 



1</>^W = {cP^(^),H} = '^(D^Es^.^iX^H)- D^,X^E^^,)n) . (6.9) 



On the other hand, a spatial differentiation of (|6.8| ) yields 



^^^r = = u;^^D^-i?B(0)(XeW) . (6.10) 



Superficially, the spatial differentiations do not commute with the time derivatives; compare ( |6.9D 
and (|10D. But we shall see that this is an illusion. First of all, it is easy to localize potential 



problems to the boundary dT,. Away from the boundary, the two expressions ( |6.S| ) and ( 6.10D 



fully agree. However, being interested in the boundary dynamics, this does not quite satisfy us. 
Let us check that they also agree on the boundary. Clearly, the expressions are singular at the 
boundary, so the only way to extract meaningful information is to prepare both the expressions 
with a general test function r]A(x). It is not difficult to see that the annihilation principle sweeps 
away any differences between the two smeared expressions: 

J^p{x)d^x r^A{x) [^,I)^]<A^(x) = . (6.11) 
6.3 Floating Poisson bracket 

The floating Poisson bracket is defined via the induced functional derivatives 



In the second equality, we used the identit y (^.17] ). Let us check that the floating Poisson bracket 
( |6.12 ) becomes the full Poisson bracket ([4.5|), by use of the annihilation principle when F and 



G are both local functionals with no external dependence. This is perhaps best seen from the 
second right hand side of ( 6.12| ). The idea is now to obtain a restrictive suitable form by recasting 



all the derivatives onto smooth objects. In the case at hand, the higher functional derivatives 
are both smooth functions of x. One may deduce that the required form can only be obtained 
for terms (M, N) when at least one of the indices M, N are 0. This reproduces precisely the full 
Poisson bracket ( ^.51 ). Technically speaking, for more general functionals F(z) and G{w), the above 
truncation takes place on parts of the functionals where both functionals are evaluated by the inner 
method. If at least one of them is evaluated by the outer method, then no truncation is carried 
out (although a truncation may take place at a later integration). This is particular the case for 
two functions F{z) and G{w). 

The advantage of the floating formulation is at least two-fold: First of all, the full Poisson bracket 
can be written in a more compact manner. The second reason is that the floating Poisson bracket 
manifestly commutes with the spatial derivatives, cf. ( ^.18 ). This is quite useful. For instance, 



consider as previously, the Hamilton equations of motions, but this time with the floating bracket 
as the Poisson bracket. We get as before 



|<^^ = {<P^,H} = '^E^^.^iX.n) . (6.13) 
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However, this time the spatial derivatives commute manifestly with the time evolution 



= {^^(^\h}-D''{cI)^,H} = 



(6.14) 



At the end of the day, it always boils down to the full Poisson bracket (4^), when all the variables 
are integrated out (and the annihilation principle applied). We voluntarily throw in a lot of formal 
zeroes in the floating Poison bracket dressed up in the above "divergent-looking" disguise. This 
gambit enables us to write the full Poisson bracket in a very compact manner. 



As perhaps the most important point, let us note that the Hamilton equations of motions (6.13) 
follows from extremizing, in the sense of ( 5.15| ), the following natural action: 



dt 



pd'^x 



^^^{x,t)LOAB<}>'^{x,t) -nix,t) 



(6.15) 



One may define, a more general floating a-bracket, for a > 1, where 

{F(.),G(u;)}(.) = j^pdxX, (-)^^A(^^ s<t>^^x) 

, i-Dr'xfjx) 5F{z) _^^ {-D)^'xfix) 6G{w) 
^ X,(x) 5,^^W(x) X,(x) 6<pm)(x)'^ ' 

In this paper, we will merely view the floating (a > l)-brackets as a curiosity, which nevertheless 
has some relevance when addressing time-order issues, see next Section |^. It coincide in the inner 
sector with the above floating (a = l)-bracket. And it satisfies the Jacobi identity. We give an 
independent proof of the important the a = 1 case in the Appendix ^ and leave the general case 
a > 1 to the reader. The a-factor slows down the e-convergence process, but viewed as an isolated 
issue, it does not jeopardize the convergence. 



7 Time Order and Quantization 

One might suspect that the boundary terms are related to the time order prescriptions. We can 
give some heuristic arguments which points in that direction. Consider first the following totally 
ordered, antisymmetric and transitive time order prescription for two spacetime points (x(i),t(i)) 
and (x(2),t(2)) in unphysical space-time IR'^ x IR. 

{X(l),t(^l)) -< {x^2),t{2)) ^ {x(^i),X(^2) G S A t(i) < t(2)) V ^ S A X(2) G 

V (2^(1), X(2) ^ S A |dE(2;(i))| > \dj:ix^2))\) ■ (7.17) 

This prescription has as a consequence, that the boundary dT, is assigned to the infinite past, so 
there effectively is no spatial boundary. (In a similar manner, one may link it with the infinite 
future.) As we shall see below, the equal-time relation ~ plays a crucial role, so let us define it 
properly: 

ix{i),'t{i)) ~ {x(2),'t{2)) ^ (a;(i),X(2) G S A = t^2)) 

V {xil),X(^2) ^ S A |dE(x(i))| = |(is(x(2))|) . (7.18) 
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We can now give the time order prescription Te for n operators = , ) , where i = 1, . . . , n. 



^(n) 



^(1) 



F, 



(^(1)) 



(7.19) 



where vr € 5^ is the unique permutation {1, . . . , n} — > {1, . . . , n}, that satisfies 

Vi = 1, . . . ,n - 1 : 7r(«+l) < 7r(i) ^ (2;(,r(^i)), i{7r{ifi))) ?^ (^^{TrW), i{7r{i))) • (7.20) 

Let us now consider an equal-time shce t = t(o)) ^-c in the traditional sense of the word "equal- 
time", as we did in the previous sections. We shall suppress the time coordinate in the following. 
We have 

rs[F(x),G(y)] = 



1e(x) lE(y) + l{o} [mx^i))\ - |o!E(x(2))|j (1 - ls(x)) (1 - lE(y))J [F{x), G{y)] 
« Mx) ls(y) [F{x),G{y)] . (7.21) 

In the wavy equality ~, we neglected a contribution from a spatial hypersurface of dimension d — 1, 
and hence of Lebesgue measure zero. Let us regularized this as 



T^[F{x),G{y)] = Xf(x) xf(x) [F{x),G{y)] , 
for some positive power /5 > 0. When time ordering the spatial derivatives we get 



Dfi^Dfy)T^[F{x),G{y)] 



Df^){x^{x)F{x)),Df'^^ (xf(y)G(y) 



(7.22) 



(7.23) 



(This equation should be understood as follows: On the left hand side the time order is every- 
where present while we form the quotient of differences for the derivative Di. In particular, it is 
present before we actually take the infinitesimal limit to produce the derivatives Di. On the right 
hand side the time order does not recognize how the spatial derivatives earlier were produced. 

only sees the result: an operator depending on one space point.) This should be compared to 
the corresponding property of the floating (a> l)-Poisson bracket 

1 1 



(a) X^i(x) X'^^y) 



{Z)5 (xri(x)F(x)) ,Df^^ {x^\y)Giy))} 



(a) 

(7.24) 

This carries some evidence, that we should translate the floating (a > l)-Poisson bracket into the 
commutator with the above perculiar time order prescription, with /3 = a — 1: 



-/3=Q-1 



[Fix),Giy)] 



{F(x),G(y)}(,) . 



(7.25) 



Although the exact value of a > 1 should be taken with a grain of salt, let us compare this behavior 
with the behavior of the floating (a = l)-Poisson bracket: 

Dff)Dg){F(x),G(y)} = {l?jy(x),Z?J)G(y)} , (7.26) 

This corresponds to the commutator with the usual time order prescription. It is also interesting 
to compare with the corresponding property of the bulk Poisson bracket, cf. ( 5.19 ), although it of 
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course does not satisfy the Jacobi identity and is therefore not expected to play any leading role at 
the level of quantization: 

So we have here presented two physically different, but both consistent, time orderings. One gov- 
erned by the floating (a> l)-Poisson bracket, but with some of the "equal-time" surfaces wrapped 
up along the boundary dT,. And another system governed by the floating (a = l)-Poisson bracket 
with the boundary being a true spatial boundary for the system. Although the above analysis 
clearly may be criticized wrt. 1) the order of various limits taken, 2) the omission of the role played 
by the annihilation principle, and 3) its disregards of further ordering issues (like *-product), it 
tends to confirm the importancy of the boundary terms of the Poisson bracket, and that they 
should not be discarded in a full treatment of a quantum field theory with a spatial boundary. 

8 Conclusions 

In this article we have 

• Reviewed the higher functional derivatives and the extended notion of differentiability of 
functionals. 

• Shown a new way to add a boundary contribution to the usual "bulk" Poisson bracket, so 
that the Jacobi identity is satisfied. 

• Given a manifest formulation of this new Poisson bracket. 

• Geometrized the Poisson bracket to an abstract world volume manifold. 

• Reviewed an embedded framework to treat Dirac delta distributions (5s(x,y) in the presence 
of a boundary 9S. 

• Introduced an annihilation principle and a floating Poisson bracket. 

• Given an action principle for Hamiltonian systems with a spatial boundaries. 

• Discussed the relation between the boundary terms in the Poisson bracket and the choice of 
time order in a heuristic manner. 

Acknowledgements. I would like to thank R. Jackiw, K. Johnson, J. Cruz, J. Pachos, M.I. Park, 
R. Schiappa and J. Song for discussions. The research is supported by the Danish Natural Science 
Research Council(SNF), grant no. 9602107. 

A Various Identities 

First of all, let us prove that the Heisenberg algebra reproduces the algebra ( |3.2D of partial and 
spatial derivatives. This follows from 

j2 ^P4(;.) d'' = Pa e'^^e^y = Pa e^^ e"^ e^^^'^^^ = e^yPA{q)e''y 

k,n>0 



23 



E E 



{qy)"" (dy) 

^ m! (n — m)l 

m>0 k,n>m ^ ' 

„k„,n mm(fc,n) 

E ^ E 



]PA{k-m)'l 



k—m 



n 



n 
m 



A{k-m) 



(A.l) 



k,n>0 ' 771=0 

Next, let us check the identity (2.15). The proof goes as follows: 



Qk 



EA{k)f 



exp 



EA{q)f S 
Jd'^y e9yEA{y)fe 



PA{q)f exp 



d- 



9=0 

A 



9=0 



Jd'^y PA{y)f e^y^-*^^ 
PA(k)f d'^d^^ 



(A.2) 



Let us note the following consequences of the Heisenberg algebra 

e'^^PAiVA) = PA{yA + y)e^y, 
EAivA) = EA{yA + y)e^y, 
EA{yA)EB{yB) = EA{yA + yB)PB{yB) = PB{-yA)EA{yA + yB) 

Similary, we mention 



EA{yA)e^^ = EAivA -y) , or equivalently EA{k)d'^ = EA{k~n) 



(A.3) 



(A.4) 



It is worth pointing out the case A; = 0, n 7^ 0, which in words says that the usual A: = Euler-Lagrange 
derivative of a total derivative term is identically zero. This is hardly surprising. Also note that 
the usual A: = Euler-Lagrange operator in this language reads 



EAik=o)f = EA{q = 0)f = J dS EA{y)f . 



(A.5) 



B Proof of the Jacobi Identity (Commutative Case) 

With the above machinary working, we can give a hopefully readable proof of the Jacobi identity. 
Consider three functions /, g and h. Using the fact that the usual Euler-Lagrange derivative cannot 
"feel" a total derivative term, we have 



{/, {9, h}} = PAia)f ^""^ 5"i^B(0){^, /i}{0) + d'EAi^o)! ^""^ PBib) P'i^c{o)5 '^^^ PDid)h 

+ Pc{c)9 d''En(o)h - {g, /i}(o)] - £^A{o)/ Eb(o){9, h}{o) 
= ri(/, g, h) + r2(/, g, h) + r3(/, g, h) - T^if, g, h) - T,{f, g, h) - {g ^ h) , (B.l) 

where we have introduced a shorthand notation for the following five terms 



Tiif,g,h) 
T2{f,g,h) 
T3{f,g,h) 

n{f,g,h) 



PB{b)Ec(o)g Efj^Q-jh 



d'^E^ 

d^EA(o)f uj^" PB{b)d'^Ec{Q)g uj^^ PD{d) > 



,AB D D „ , fD Qc J 
,AB D fld 17 „ , .CD 



^A(O)/ PB(b)Pc{c)9 "^""^ d'^Ej^f^Q-^h 



d^EA{f))f io^^ PB(b) Ec(o)9 "^^^ E 



CD 



D(0) 



EA{o)f ^^"^ i-d)' [PB(b)Ecio)9 ^""'^ Er)^o)h ■ 
The Jacobi identity, containing 30 Tj-terms, now follows from the fact that 

T2{f,g,h) = T2{h,g,f), Ti{f,g,h) = T3{h,gJ), T^{f,g,h) = n{h,g,f). 



(B.2) 
(B.3) 
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The first equation is trivial and the next two equations follows by rewriting in terms of Fourier 
transforms 



T3if,9,h) 



T,{f,g,h) = PA{qA)f^^'' exp 



d 



d 



dqA dqs 



q=0 



d'^'y PA{yA)f t^^^ e^^y^-y^^ 



\PB{<lB)Ec{<lc)g ET){(lT))h 

PB{yB)Ec{yc)9 ED{yD)h 
e-^y^PB{yB+yc)Pc{yc)9 '^^'^ ED{yD)h 



I d^^y PA{yA)f 

J d'^'y PA{yA)f o;^^ e^^y^-y^^PB{yB)Pc{yc)9 a;^"" e^^y^+y<^-y'^^ED{yD)h 



cxp 



d 



d 
dqB 



d 



d 
dqo 



Ec{qc)9 ^'^'^ PD{qD)h 



Ti{f,9,h) = exp 



EAiqA)f oj^"" PB{qB) exp 
j d^^y e^y^EA{yA)f u^^"" PBiys) e^y'^Ec{yc)9 ^^"^ PD{yD)h 
I d'^'y e^y^EA{yA)f o;^^ e^^y'^-y^^PB{yB+yc-yD)Pc{yc)9 o;^"" PD{yD)h 

J d^'^y e^^y^+y^-y^^EAiyA)/ w^"" e^^y'^-y^^PB{yB)Pc{yc)9 ^""^ PD{yD)h , 

exp EA{qA)f OJ^"" PB{qB) Ec{qc)9 ^^"^ ED{qD)h 

L dqB\ 5=0 

= j d^S e^y^EA{yA)f o;^^ PB{yB) Ec{yc)9 En{yD)h 

= J d'^y e^y^EA{yA)f u;^^ e-^y^PsiyB + yc) Pc{yc)9 En{yD)h 

= J d'^'y e^^y^-y^^EAiyA)/ o;^^ e-^2/cp^(y^) Pc{yc)9 o;^^ ED{yD)h , 



n{f,9,h) = EA{qA)fu;^'' exp 



-d 



dqB 

dys 



g=0 



PB{qB)Ec{qc)9 ^'^^ ED{qD)h 

PB{yB)Ec{yc)9 ED{yD)h 
e-^y^PB{yB + yc)Pc{yc)9 ^^"^ ED{yD)h 



d'^^y EA{qA)f cu^^ e 
= j d^S EA{qA)f o;^^ e 
= I d^^y EA{qA)f oj^^ e-^y^PB{yB)Pc{yc)9 ^""-^ e^^y<^-y^^ ED{yD)h 
We have used the following shorthand notation for the integration measure 

d^'^y = d%A d'^yB d%c d'^yo , 
and we have performed the following change of integration variables 



(B.4) 



(B.5) 



(v'a] 




( 1 












(yA \ 


v'b 







1 


1 


* 




yB 


y'c 










1 







yc 


\ v'd ) 












1 J 




[yo J 



(B.6) 



which has Jacobian equal to 1. 



□ 
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C Calculus of Words 



To be self-contained, we will here give a short treatment of the calculus with ordered index-structure, 
merely giving the main definitions and formulas. 

C.l Words 

An ordered tuple K (or a positive word) takes the form 

K = (A:i,...,fc|^|)G {!,..., . (C.l) 

Here we have employed a d-letter alphabet {1, . . . ,d} and \K\ € INq denotes the length of K. The 
transposed word is = (fc|/^|, . . . ,ki). We define the (non-commutative, associative) sum of two 
tuples K = {ki, . . . , A;|^|) and L = {ii, . . . ,i\L\) the concatenation 

K + L = {ki,...,k\Ki,h,...,iiL\) ■ (C.2) 

Obviously, the empty index set is the neutral element. Formally, we can define negative words 
—K = —{k^K\) — ■ ■ ■ — (ki). Moreover —{K + L) = —L — K. We define a left subtraction —K + L 
as the unique solution lo K + {—K + L) = L. If the left subtraction —K -|- L is a positive word, 
we say that K < L. Clearly, {K + L)* = L* + KK 



There exists another partial order between two arbitrary positive words K and L. Namely, define 
that K preceed (or is equal to) L, written K ^ L, ii we can obtain K from L by deleting some 
(possible no or all) elements in L. Said in a mathematical precise manner, there exists a (strongly) 
increasing index function vr : — > called an (orderpreserving) embedding, 

such that 

h = (^Tr(l) , ■■■ , k\K\ = K{\K\) ) 

7r(l) < ... < 7t{\K\) . (C.3) 

In the affirmative case, we define the subtraction L ^ K as, the tuple of deleted entries. More 
precisely, in this case there exists a unique (strongly) increasing index function, called the comple- 
mentary embedding, 

TT^ :{!,..., , (C.4) 

such that the images of tt and vr'^ are disjoint, i.e. 

^({i,...,|i^|}) n ^^({i,...,|L|-|if|}) = 0, 

^"(1) < ... < 7r^(|L|-|i^|) . (C.5) 

Then the subtraction L ^ K is defined as 

L^K = (4c(i), . . . ,£^c(|i|„l^l)) . (C.6) 

We stress that the embedding vr is not necessary unique. Therefore L ^ K depends on the em- 
bedding vr. For instance, in the entanglement formula (4.7), one should sum over all possible 
embeddings. A closed expression for the degeneracy d{K ^ L) of imbeddings vr is not known to the 

author. By definition d{^^L) = 1. Note that K ^ L implies L -\r K ^ L. In the affirmative case 
the pair K and L ^ K is called [^] an unshujfle of L. There exist 2l^l unshuffles of L. Furthermore, 

(L H- KY = L*^K\ 
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K < L < N =^ N^L = {N^K)^{L^K) . (C.7) 

A shuffte K^M of two positive words K and M is defined as the opposite of an unshuffle in the 
sense that it is a solution X to X ^ K = M. Clearly, the number of shuffles for fixed K and M is 

^ ^^^Im]"^^^ ^ ' '^^^ number of solutions X to L -\r X = M for fixed positive M, L with M ^ L is 
d{M<L) . 

C.2 Alphabets of Operators 

A d- letter alphabet of operators (or more generally, of associative abstract algebra elements), is just 
d operators A = {Ai, . . . , A^). The sum and the product {i.e. usually the operator composition) of 
two alphabets are defined letterwise 

A + B = {Ai + Bi,...,Ad + Bd) , AB = {AiBi, . . . ,AdBd) , (C.8) 

respectively. 

C.3 Words of Operators 

If we have an alphabet of operators A = {Ai, . . . , A^) we can form words of operators 

= Ak,...Ak^^^ . (C.9) 

We invoke the convention that Ak = A*-^'^ denotes the transposed word. The empty word operator 
A'^ = 1 is the identity. Operators for non-positive words, i.e. for words containing negative letters, 
are declared to be zero. Concatenation leads to a non-commutative product 

between words (and A^Am = Am+n for the transposed). It coincides with the letterwise multi- 
plication. But this is not the only associative product of words. Shuffling leads to a commutative 
*-product 

where the sum is over possible shuffles K^L. The definition of the *-product is extended by W- 
bilinarity. The concatenation product and the *-product coincide for commutative alphabets. The 
following binomial relation is a consequence of the special features of unshuffles, cf. (|C.8| ) and ( C.9 ), 



[Ai,Bj] =0 {A + B)^ = J2 A^'^B^^^'^ . (C.12) 

M=0 

We can also define a non-commutative, associative sum by the following binomial expression 

This sum and the usual sum do not coincide in general, not even for commutative alpha- 
bets. But see the equation ( p. 22 ) below for further relationship. 
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Note that the definition ( C.ll| ) makes no use of the algebra multiplication. It only needs a vector 
space with a basis of vectors indexed by words K. The definition of the #-sum (|Cl3|) needs 
an algebra of words , but it is irrelevant whether A^ is a composite object of more elementary 
letters or not. The same remark could be made about formulas ( ClOj ) and ( p.l2| ), that in the 
minimalistic interpretation becomes definitions. Even in the case where there exists a letterwise 
algebra multiplication, we will often use the binomial formula ( |C.12D , also known as a convolution, 
which makes sense even for mutually non-commuting alphabets. 

We have that -{A#B) = and 

Sometimes we will also need to define = J2k\^^ ^\k\i |TjT' 

C.4 Functions of Operators 

For an analytic function /(x) = Z^J^Lo^n^"") define f{A) = J21^=ii)(^\N\A^ ■ In particular, the 
exponential of an alphabet is 

Af=0 I I' 



We have that 



[Ai,Aj] =0 exp(^) = exp(Ai)...exp(^d) , 

[Ai,Bj] =0 ex.p{A + B) = exp(A) exp(5) . (C.16) 

Also we have the important orthogonality relation 

A/=0 

This property leads to the vital relation e^e^^ = 1 even for a non-commutative alphabet. Similary, 
we have ((-^)*#^)^ = 0^. 

In practice, we only use f{AB) for a product of two alphabets. (Dummy indices usually come in 
pairs.) It is convenient to define a "dot product" notation 

oo 

fiA-B) = ^ a|^|yl^S(^') , (C.18) 

implementing a transposition of one of the alphabets. Also /((— • B) = f((—A)B). We have 
the following inversion relation for the concatenation product 

[Ai,Bj] = exp[{-A) ■ B]exp{A- B) = 1. (C.19) 

This should be compared with the inversion relation for the the mixed case of an implicitly written 
concatenation product and a *-product 

[Ai,Bj] = exp\{-AY-B] *b exp{A ■ B) = 1. (C.20) 
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We have the following distributive laws 



[Ai, Cj] 
[Bi,Cj] 









f exp [(-A)* • 

exp[^-(B#C)] 
exp [{A + B)-C] 



exp((- 


Af ■ B) exp((- 


Af 


exp(A 


■C) exp(^-5) 




exp(^ 


• C) *c exp(5 


■C) 


exp(S 


• C) *c exp(^ 


■C) 


#" and the usual sum 


U j » 



loosely speaking in average. More precisely, 



[Ai , Aj 







exp[^-(S#C)] = exp[^-(B + C)] 



(C.21) 



(C.22) 



C.5 Trace and Fourier Analysis 



We define a trace on the vector space of words, which is constructed from two mutually commuting 
freely generated alphabets A and B: 



J d'^A d^B e(-^)*-^ *B A^ Bm = Ti{A^ Bm) = \N\\ 5^ , 



(C.23) 



and extend by (T-bilinarity. As the first equality suggests, we will sometimes use a suggestive 
notation for the trace borrowed from the Fourier analysis in the usual commutative case. One can 
take this analogy quite far. We do not give any meaning to the position of the measure d'^A d'^B, 
i.e. it is taken to commute with everything. A theoretically perhaps more convenient form is 



J d'^A d'^B exp [{-A - Aj ■ {B#B') 



= 1 . 



(C.24) 



Prom here it follows trivially that the integration measure d'^A d'^B is translation invariant under 
A^A + A',B^ B#B' (but not under B B'#B\). 



C.6 Differentiation of Words 

d _ / a 



Consider the differential alphabet ^ = (g^, . . . , g^^) of freely generated associative algebra ele- 
ments Ai. Let us define differentiation at the level of words {i.e. not letterwise), as 



d 



iKWdA^ 



dA^ 



[1] 



(C.25) 



where the sum is over possible embeddings tt. Extend the definition by (T-bilinarity. One of the 
main motivations behind this definition is to implement the Taylor formula 



[Ai,Bj 







exp 



B ■ 



d_ 

'dA 



f{A) = f{B + A) 



The composition of derivatives is described by the *-product 

1 d 1 d ^ 1 d 1 



iKWdA^ iLWdA^ 



E 



d 



d 



\N\\dA^ ~ \K\\dA^ \L\\dAL 



(C.26) 



(C.27) 



N= K#L 

As consequences, the derivatives are associative and commutative wrt. composition. They enjoy 
the following properties 

d 



[Ai,Bj 







exp 



i-Bf 



d_ 

dA 



exp 



B 



dA 



1 , 
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[Ai, Cj] 
[-Si, Cj] 








exp 



d ■ 








C 


d - 


91. 


= exp 


a4. 


exp 


'dA 



. (C.28) 



We can of course implement the concatenation product for the differential. However, in practical 



calculations it plays no role. The derivatives ( C.25 ) do not satisfy the Leibnitz' rule. 



The are other kinds of differential alphabets. In our case, we have the covariant derivatives Di that 
act on words {—Y)k according to 



[Di,{-Y)K] 



\K\ (-y)_(,)+i, if (i) < K 
otherwise . 



(C.29) 



(One can consider the words {—Y)k as originating from an associative algebra alphabet {—Y^, . . . , —Y'^ 
that behaves non-associatively wrt. the DiS, but it is unnecessary.) The covariant derivatives sat- 
isfies Leibnitz' rule on functions: 



(C.30) 



M=0 



This can be recasted into the Taylor-like form e^'y{fg) = e^'^ f *y e^'^g. 



D Proof of the Jacobi Identity (Non-Commutative Case) 

We have 



{f,{9,h}} 



+ Pc{C)9 ^ 



CD r^C, 



B{B) 
AB 



D ED{o)h - {g,h}^Q) - EA(Q)f u; EB(o){g,h}(^Q^ 



= Ti{f,g,h)+T2{f,g,h)+T3{f,g,h)-T4f,g,h)-n{f,g,h) - {g ^ h) 
where we have introduced a shorthand notation for the following five terms 



D{D)' 



(D.l) 



,CD 



^D(O) 



T2{f,g, h) = D"EA(^-)f PB{B)Pc{c)g ^'^'^ D^Enio)h 
T3if,g,h) 



D^^Ea^o)/ '^"^ PBiB)D''Ecio)9 



,AB 



CD 



n{f,g,h) ^ D^E^^o)! PBiB) Ecio)9u;''''' E_ 



AB 



.CD 



D(D) 



n{f,g,h) = Ea^o)/^^'' {-D 



D(0) 

PB{B)Ec{o)g ^'^^ E 



D(0) 



(D.2) 



Here we have chosen to use the same index symbol A, B, C and D to label the indices of the fields 
(p and the words. It should not lead to any ambiguities, and it hopefully becomes easier to grasp 
the index structure. The Jacobi identity, containing 30 Tj-terms, now follows from the fact that 

T2{f,g,h) = T2{h,gJ), Ti{f,g,h) = n{h,g,f), mf,g,h) = n{h,g,f). (D.3) 

The first equation is trivial and the next two equations follows by rewriting in terms of Fourier 
transforms 



Tiif,g,h) 



PaJ u:^"" exp 



D 



_d_ 

dqA 



exp 



d 



dqi 



PsEcg oj^^ Enh 



q=0 
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T3{f,9,h) 



n{f,9,h) 



Ttdcba PaJ t^^"" e^-J^-^e^-^)'-^^ [PBe^~''^' y<^ Peg t^^^ EdH 

Ttdcba PaI u;^^ e''-y^e(-''^' y^ [e^-'^-'^^^'-yo P^Pcg u;^'' Eoh 
TvDCB'A PaJ c^^"" e''-y-^e'--''^' y'BPBPcg cc;^^ e^'^^e^-^^^-J'se^-f^i^i^/i 



exp 



D ■ 



_d 
dq 



B 



EaI oj^^ Pb exp 



D ■ 



_d 

dq 



D 



Ecg Poh 



q=0 



TvncA e'^-y^EAf o;^^ Tr^ Pse'^-y^ e^-'^^'-y^ Peg o;^^ Poh 
TvDCA e^-y^EAf uJ^^ Trs e^^^^'^^yvo fA-D-quf-yc p^p^g ^CD 

TrccA e^-y^e^-^^'-y^^e^-y'BEAf Tr^, e^y° e^-^^'-y^ PsPcg oJ^^ Poh , 
^ ^ EaI uj^"" PBEcg i^^"" EdH 



exp 



D ■ 



dqB 



q=0 



TvDCA e^'-y^EAf o;^^ Tr^ Pse^-^^^-^^PcS ^^"^ Eoh 
TTDCAC^'-y^EAf u;^^ Ttb e^~''-''^^'-y^ PePcg u^""" E^h 



Tr^cA e^-^'^'-y'='e^-y'BEAf uj^^ Ttb' e^'^'^' y^ PePcg ^^'^ Eoh 

^ ^ \PBEcg io^^ Eoh 



EAfu;^"" exp 



i-ny 



dqB 



g=0 



Ttdcba EaI e^-^fy^ Pse^-^'^'-y^ Peg u;^^ Eoh 



Ttdcba EaJ e(-^)*-^- 



■y^PsPcg u;^^ Eoh 



Ttdcb'A EaI u;^^ e^-'^'^'-y'BPBPcg u;^^ e^-^^'-^se^-^^c^^^/, . 



(D.4) 



We have suppressed *-products among the y- variables, and shorten Pa = PA{qA)^ Pb = Pb{qb) 
. . ., Ea = EA{qA), etc. The trace can be written more suggestively as 



Ttdcba 



We have performed the following type of translation of the integration variables 



v'b 



yB#yc or y's = VB^i-yof^yc ■ 



(D.5) 



(D.6) 



Note that after the shift of integration variables the y-alphabets do no longer mutually commute. 
However, one may convince oneself that the integrations can be unwind, and we can consistently 
declare them to mutually commute also in the new variables. Finally to prove (D.3), one should 
relabel dummy variables ABCD DCBA. 



□ 



E Proof of the Jacobi Identity (Floating Type) 



We now turn to the proof of the Jacobi identity for the floating Poisson bracket, cf. Eq. ( |6.12| ). Con- 
sider the local functionals of Subsection We assume that Dip = 0. Suppressing the integrations, 
we have 



E 



A{0) 



(^^/) , ,AB p 



B{0) 



Ec{o){^eg) >^ 



CD 



E 



D 
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T{f,9,h) - (g^h), 



(E.l) 



where 



u^^ exp 



PB(B)Ec(Q){^eg) ^ 



CD ^D{0)i'^eh) 



d 
Qqb 



PBEciX,g) u 



CD EDiXeh) 



q=0 



Ttdca e^^-^y-yce^^-y'B^^^ 

Xf 



PBe(-^y-yoPc{X,g) 



X, 



^i-D-qsY-yc 



PBPciXeg) OJ 



CD EoiXeh) 



X, 



e(-^y-ycp^Pc{X,g) u;^^ 



X, 



= Ttdcba e^(-^)*-^^e^^-^«:^^^ cv^^ e5(-^)*-^-e^(-^)*-^cp^p^(X^^) ^CD 

Xe 

,^-i-Dr.y,^lD.yaM^ ^ T{h,g,f). (E.2) 



X 62^ 



X, 



The ~ indicates that the equahty holds up to total derivative terms. They are unphysical terms 
living far away from the bounded physical region S, and therefore vanishing. In the last step 
we substituted y'^ = hbH^IIC- The annihilation principle will not change the fact, that the Jacobi 
identity is fulfilled, because all annihilated terms appear in pairs with opposite sign. 



□ 



Let us now turn to more general functionals F{u), G{v) and H{'w). We have 

6 



{F{u),{G{v),H{w)}} = Jp{x)d''xp{y)d^yF''{a 



= I p{x)d''x p{y)d''y TiF, G, H) - {G{v) ^ H{w)) . 

Here we have applied the following shorthand notation 

5[\.(.-r)F(»)] _ ^CD s[X,Xy)II{w)\ 



F^(x) 



6^^{x) Xe{x) ' 



and 



T{F,G,H) ^ F^{x){-D^,.)^ 



= F'^{x){-D^,)Y 



X,{x) 



X,{x) 



X,(y) 64>^{y) 

5[X,{y)G{v)\ 



(E.3) 



(E.4) 



d^B(B){x) 5(l>G{y) 



d 



X 



Ti{F,G,H)+T2{F,G,H) , 



p{y) ^^(^) 



G{v) 



(E.5) 



where the last equality will be explained below. We distinguish between the so-called inner j- 
terms j = 1, . . . , s, where the spatial i?^^) -differentiation are applied on the Pc'-derivatives of the 
G-functional before the partial derivative Pb, and on the other hand the so-called outer j-terms. 
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where the order is the opposite. For each inner j-term, we may write G{v) = J p{v(^j-^)d'^V(^j^ djiv)- 
Together, the diagonal piece of inner j'-terms becomes 



T,{F,G,H) ^ 



-^W) ^(^^ Pb(b) (-^(-u))) Pc{c)9m\,^^^^^ 



1y /'(^) 



p{y) 

The y-integration may be exphcitly performed in the diagonal Ti-term: 



(E.6) 



(E.7) 

It may now be treated similarly to the local case discussed in equation ( [E.2D . The rest of the terms 
appearing in ( |E.5| ) can be organized so that they are manifestly symmetric in F and H, and hence 
do not effectively contribute to the Jacobi identity: 



T2iF,G,H) ^ F'^ix) 



p{x) p{y) 'B{B)^C(C) 



H'^iy) ■ 

(E.8) 

(The prime ' indicates that the above inner diagonal j-terms should not be included in the T2-sum.) 

□ 



F Realization of Derivatives 

In this section we will like to incode in an alternative manner the information about which part 
of an expression that are hit by a derivative. Instead of the usual practice of indicating this with 
arbitrarily many arrows, we define a linear chain of operators containing the same information. 
Consider some expression, where upon the derivatives Di act. Notation: For simplicity, we assume 
that the covariant derivatives commute, i.e. that the curvature vanishes. First of all, we assign 
names a to the derivatives Di -D(o)i so that we may distinguish them. For all derivatives 
let us introduce a scalar elements a(^a) and a vector element 6*^^^ that obey 



1(a) 0(0) - 1 , ^(a)^(") ~ ^ ' ^(a)"(a) 



(F.9) 



The notions scalar and vector refer to their properties under coordinate transformations. We have 
the following commutator relations 



a(a))«(/3) 



0, 



^{°)'^(/3) 



0, 



«(a)'^(a) 







"■(a)' ^(/3) 



25^p IP^^^ . (F.IO) 



Furthermore, we introduce a Berezin-like integration J db(^a)i- I^i the following, / and /' denote 
functions which do not depend on a(a) and but possibly on other a's and 6's. We will assign 



the value 



J db(ayi h'^^^ f = §j f , J dh^a)i f 



. 



(F.ll) 
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We emphasize that we have not defined the integral / db(^a)i if there are a(a) elements present in 
the integrand. Let us also declare 



f 







, 



f, /' 







(F.12) 



Now we are ready to define the transition from derivatives to a's and 6's. The derivative D(^a)if 
may now be written as 



This on the other hand is evaluated by declaring 



db 



{a)i 



He) 



(F.13) 



(F.14) 



In this way we may imagine an expression with an arbitrary even (non-negative) number of a(^a) 's- 
For instance, one may check that 



db. 



{a)i 



f ata) 9 /' a(Q) a' aia) f" = ff'f" Di{gg') 



(F.15) 



Here we see that the a's plays the role of "start and stop signs" for the action of the derivatives. 
The Leibnitz rule is implemented via the assertion ( |F.14| ) . An integration by part of an expression 
may be performed by cyclicly moving an a(^a) from one end of the expression to the other end 
under an additional change of sign. Let now s and s' be two singular expression, cf. point Bl in 



Subsection 5.5. Some of the consequences of the annihilation principle may now be stated as 



1(a)) 







(x(c,)) a(c,), 5{x(^a)-y) s'{y) 



. 



(F.16) 
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